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We consider the excitation of the nuclear transition 2299 Th - 229mr ph near 7.6 eV in singly ionized 
thorium via an electronic bridge process. The process relies on the excitation of the electron shell 
by two laser photons whose sum frequency is equal to the nuclear transition frequency. This scheme 
allows to determine the nuclear transition frequency with high accuracy. Based on calculations 
of the electronic level structure of Th + which combine the configuration-interaction method and 
many-body perturbation theory, we estimate that a nuclear excitation rate in the range of 10 s~ 1 
can be obtained using conventional laser sources. 
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The nuclear transition between the low-energy isomeric 
state and the ground state in the 229 Th nucleus is very 
interesting due to the possibility to build a very precise 
nuclear clock and its high sensitivity to a hypothetical 
temporal variation of the fundamental constants 0, 
The value of the transition frequency is known from 7- 
spectroscopy with an uncertainty that is many or- 
ders of magnitude higher than the presumed linewidth. 
The presently most precise value of the transition energy 
is 7.6(5) eV @, placing the transition in the vacuum- 
ultraviolet spectral range. Further investigations are 
therefore needed in order to allow a direct optical excita- 
tion and detection of the line and to gradually increase 
the spectroscopic resolution to the level adequate for an 
optical clock [lj. Since the energy of this nuclear transi- 
tion is in the same range as resonances of the outer-shell 
electrons in thorium ions and chemical compounds of Th, 
the electronic environment may have a significant influ- 
ence on the nuclear transition rate In this work we 
suggest a new experiment to study laser excitation of the 
nuclear 229 ffTh - 229m Th transition, making use of the 
dense electronic level structure of Th + . 

In a recent article we have studied the effect of 
atomic electrons on the nuclear transition from the iso- 
meric 229m Th state to the ground 229 3Th state in 229 Th+ 
due to an electronic bridge (EB) process. Here we con- 
sider the process of the nuclear g — >■ m excitation via an 
EB process. Atomic units (K = \e\ = m e = 1) are used 
unless noted otherwise. 

Experimental scheme. The EB process considered 
in this work can be represented by the Feynman diagram 
in Fig.[T] In the following we assume a resonant character 
of this EB process. For this reason we take into account 
only one Feynman diagram given by Fig. [1] which mainly 
contributes to the probability of the process. Other di- 
agrams which can be obtained from this one by permu- 
tations of the photon lines will be neglected. We sup- 
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FIG. 1: Feynman diagram of the considered two-photon elec- 
tronic bridge process. The single and double solid lines relate 
to the electronic and the nuclear states, respectively. The 
dashed lines are the photon lines. 



pose that the initial state i is the electronic ground state 
and the initial and final electronic states are the same. 
Further we assume a two-photon excitation of high-lying 
states k. For the first excitation at u>i, it seems con- 
venient to use a continuous laser tuned to the transi- 
tion from the Th + ground state (6d 2 7s)J = 3/2 to the 
(6d7s7p) J = 5/2 state at 24874 cm" 1 [8]. This electric 
dipole (El) transition is known as the strongest emis- 
sion line in the Th + spectrum and it can be excited ef- 
ficiently in collisionally cooled ion clouds in a radiofrc- 
quency trap @. As a result, only the (6d7s7p) J = 5/2 
state needs to be accounted for in the sum over interme- 
diate states n in Fig. 1. We also note that the moni- 
toring of the fluorescence resulting from excitation of the 
(6d 2 7s) J = 3/2 -4 (6d7s7p) J = 5/2 transition allows us 
to detect the excitation of 229 Th + ions from the nuclear 
ground state to the isomeric state, making use of the dif- 
ferences between the hyperfine structure characteristics 
of electronic transitions in 229m Th + and in 229 fTh + [lj]. 

The second excitation step with frequency uj gives rise 
to the excitation of higher- lying states k. Their decay to 
the ground state is accompanied by the nuclear g — > m 
transition and leads to the appearance of isomeric Th + 
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ions in the electronic ground state. For this excitation 
scheme, energy conservation implies that u)\ + u> = cjjv, 
where w\ and w are the frequencies of the incident pho- 
tons, and ljm = E m — E g is the nuclear transition fre- 
quency as determined by the difference between the iso- 
meric nuclear energy E m and the ground-state energy 
E g . 

This excitation scheme offers the prospect to determine 
the nuclear transition frequency lom with the accuracy 
afforded by high-resolution laser spectroscopy. In an ex- 
periment based on this scheme, one would use a widely 
tunable laser source to produce narrow-bandwidth radi- 
ation with variable frequency u>. If the laser frequency 
is scanned over the nuclear resonance, the probability 
of excitation to the isomeric state exhibits a resonance 
peak. In an ion-trap experiment with collisionally cooled 
229 Th + ions, the width of the resonance is determined 
essentially by Doppler broadening and by the combined 
linewidth of the employed laser sources. If one assumes 
saturated excitation to the state n and a fixed detuning of 
oj relative to the electronic transitions n —> fc, the prob- 
ability of the EB process shown in Fig. 1 is proportional 
to the spectral intensity of the laser field at u. 

Based on the result of Ref. Q we assume that the most 
interesting range of the electronic excitations is around 
7.6 eV w 61300 cm -1 . Following this assumption, we 
expect that in the sum over the intermediate states k 
the atomic energy levels lying between 60 000 and 64 000 
cm -1 will give the predominant contribution to the prob- 
ability of the EB process. Unfortunately, these energy 
levels are not yet identified experimentally. Therefore all 
the following results are based on ab initio calculations. 
As we noted in |7J, the achieved accuracy of calculations 
of the high-lying states of Th + is at the level of sev- 
eral percent. This is not sufficient to reliably predict 
the resonance enhancement occuring in the EB process. 
Hence, the experimental identification of the energy lev- 
els should be considered as the next step towards realizing 
the method considered here. 

Because at present very accurate calculations are not 
needed we will make one more assumption simplifying 
the calculations. We assume that only one intermedi- 
ate state n, the state at 24874 cm" 1 , contributes to the 
probability of the EB process and that 100% population 
of this state can be achieved. As a result, the frequency 
uji is assumed to be fixed at 24874 cm -1 and the process 
which we discuss in the following can be described by the 
diagram represented by Fig. [2] where the state at 24874 
cm -1 is denoted as t and is considered as the initial state. 

Calculation. Figure [2] shows an EB process that re- 
lies on the absorption of an incident photon. As we men- 
tioned in Ref. [10] an EB process of this type can be 
effectively treated as a "generalized" electric dipole tran- 
sition from the initial to the final state. If the incident 
radiation with spectral intensity I u is isotropic and un- 
polarized, the relation between the probabilities W a b of 
a spontaneous transition a — » b and of the corre- 
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FIG. 2: Simplified Feynman diagram of the considered elec- 
tronic bridge process. Here it is assumed that the first excita- 
tion step leads to efficient population of the state (Qd7s7p) J = 
5/2 at 24874 cm" 1 which is denoted as t. 



sponding induced transition b — > a is given by [Il| 



wr a = w ab 



(i) 



The spontaneous EB process can be formally described 
by the mirror-image of Fig. [2] with an outgoing photon 
arrow. The general formalism was discussed in detail 
in two recent articles @, [l(J. Assuming the resonance 
character of the EB process the expression for Teb can 
be written as 0: 
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where A4\ is the magnetic dipole nuclear moment and 
\I g ) = 5/2+ [633] and \I m ) = 3/2+ [631] are the ground 
and the isomeric nuclear state, respectively, given in their 
Nilsson classification. Ji is the electron total angular 
momentum of the initial state and u> = u>n — u>i is the 
frequency of an absorbed photon. The explicit expression 
for the coefficient Gi was derived in [jjl and reads as 
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where uj S i = e s — £j, and e; denotes the energy of level I, 7i 
is the electronic magnetic-dipole hyperfine coupling op- 
erator, D is the electric dipole moment operator, and 7 S 
encapsulates all other electronic quantum numbers. The 
explicit expressions for the matrix elements of the opera- 
tors Ti and D are given in Ref. [T(| • The possible values 
of the total angular momentum Jk are determined by the 
selection rules of the operators 7i and D. In particular, 
in our case Ji = 3/2 and Jt — 5/2 and, correspondingly, 
Jk can be equal to 3/2 and 5/2. 

In the following we will use the dimensionless quantity 
Pmi introduced in Ref. [l(| and defined as the ratio of 
the probability of the spontaneous EB process Teb to 
the probability of the spontaneous Ml radiative nuclear 
to — > g transition T^: 
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We use the method of calculation described in de- 
tail in Ref. [7]. We consider Th + as the atom with 
three valence electrons above the closed-shell core [Is 2 , 
... ,6p 6 ] and employ the CI+MBPT approach combining 
the configuration-interaction (CI) method in the valence 
space with many-body perturbation theory (MBPT) for 
core polarization effects 12]. At the first stage we solved 
Dirac-Hartree-Fock (DHF) equations [l3| in V N ~ 3 ap- 
proximation and then we determined the 5/, 6d, 7p, 7s, 
and 8s orbitals from the frozen-core DHF equations. The 
virtual orbitals were determined with the help of a recur- 
rent procedure [lij]. The one-electron basis set included 
ls-18s, 2p-17p, 3d-16d, and 4/-15/ orbitals on the CI 
stage. 

We formed the configuration spaces allowing single, 
double, and triple excitations from the 6d 2 7s configura- 
tion (for the even states) and from the 5/7s 2 configura- 
tion (for the odd states) to the 7s-13s, 7p-12p, 6d-lld, 
and 5/-10/ shells. An inclusion of all possible (up to 
triple) excitations allows us to take into account most 
completely the configuration interaction for all consid- 
ered states. The energies and the wave functions are 
determined from the eigenvalue equation in the model 
space of the valence electrons 

H eS (E p )\<S> p } = E p \<Z> p ) ) (5) 

where the effective Hamiltonian is denned as 



Using Eq. ([3]) we obtain 



H eS (E) = H FC + 
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Here Hfc is the relativistic three-electron Hamiltonian in 
the frozen core approximation and £(£?) is the energy- 
dependent core-polarization correction. 

Together with the effective Hamiltonian H c ff we in- 
troduce an effective electric-dipole operator D e g and an 
operator (7I) c ff acting in the model space of valence elec- 
trons. These operators were obtained within the rela- 
tivistic random-phase approximation (RPA) [III, [l|| . On 
the stage of solving the RPA equations and calculating di- 
agrams for effective Hamiltonian and effective operators 
D and 71 we used a more extended basis set. It consisted 
of ls-22s, 2p-22p, 3d-22d, 4/-22/, and 5g~16g orbitals. 

Results and discussion. In the following, we assume 
that the value of ujn is between 60 000 and 64 000 cm -1 , 
which corresponds to the range 7.4-7.9 eV and suppose 
that the main contribution to Gi (see Eq. ([3])) comes from 
intermediate states lying in this range. As noted above, 
the odd-parity state (6d7s7p) J = 5/2 at 24874 cm -1 is 
assumed as the initial state t. In [1 
the largest value of the coefficient 
initial state and the intermediate state k (see Fig. 2) are 
connected by an El transition. Thus, the states t and 
k should be of opposite parity and we have to consider 
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where the quantity R s .j k is determined as 

i? SjJfc ee |<6d 2 7s, J = 3/2||7i||7aJfc> 

x { ls J k \\D\\6d7s7p, J -5/2) | 2 . (8) 

As follows from Eq. (j4]) the equation for (3 mi reads as 

Pmx « ^ f 1 - — ) G 2 « 0.012 x G 2 , (9) 
18 V J 

where we took into account that the quantity uji/ujn — 
0.4 if u N is between 60000 and 64000 cm" 1 . 

High-lying energy levels of Th + in the range from 
60000 to 64000 cm -1 were determined theoretically in 0]. 
In Table [fl we list the values of the coefficients Rs,j k 
found for the even-parity states in the frame of the 
CI+MBPT+RPA approximation. The good correspon- 
dence between experimental and theoretical levels at 
lower energies makes us confident that the level struc- 
ture is complete for the considered electron configura- 
tions. We expect that there will be more levels from 
other configurations but that the coefficients R for con- 
figurations with multiple electron excitations are smaller 
than the dominant ones calculated here. 

Using Eq. 0, the tabulated coefficients R s ,j k , and en- 
ergy values we can find Gi and /3mi for a given value of 
uin- The values of G2 and /3mi depend critically on the 
position of ojn relative to the electronic levels. This is 
indicated in Fig. [3] which shows the variation of G2 with 
ujn in the considered wavenumber range. As seen from 
this figure, G2 > 800 if a; at lies within the range of the 
electronic energy levels listed in Table HI For ujn = 7.6 
eV the (5f6d7p) J = 5/2 state at the calculated value 
60462 cm -1 yields the main contribution to the EB pro- 
cess and one obtains G2 ~ 2700 and @mi ~ 30. When 
G2 can be calculated on the basis of spectroscopically de- 
termined energy levels, it is likely that the dependence of 
G2 on wat will look somewhat different from Fig. [3] Nev- 
ertheless, we expect that the EB excitation rate depends 
on a small number of dominant channels or even only a 
single one and that the typical minimum value of the co- 
efficient G2 will remain at the same order of magnitude 
(~10 3 ). 

The excitation probability to the isomeric state due 
to resonant laser excitation with ui\ + ui = ujn can be 
estimated as follows. We assume that radiation at to is 
produced by a laser source which emits pulses with 10 mJ 
energy and with a spectral width of Aw = 2tt x 3 GHz 
at a repetition rate of 30 Hz. These characteristics can 
be achieved, for example, with commercially available 
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FIG. 3: Dependence of the coefficient G2 (see Eq. (J7J)) on the 
nuclear transition frequency ujn as calculated on the basis of 
the data listed in Table [T] 

TABLE I: Calculated energy levels (see text) in the range 
from 60000 to 64000 cm -1 and coefficients Rk,j k - is the 
difference between the energies of the excited state and the 
ground state. The notation y[x] means y x lO^. 
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frequency-doubled dye lasers as they are used in many 
laboratory applications. With focussing to a beam cross 
section of 0.1 xO.l mm 2 , the resulting time-averaged spec- 
tral intensity is I « 1.5 x 10~ 3 (W/m 2 )s. Taking into 
account Eqs. (Q} and Q and assuming /3mi ~ 30 and 
the value Tn — 6.6 x 10~ 4 s _1 [l(|, we find an excitation 
rate in the range of ~ 10 s" 1 . Significantly larger 
excitation rates can be expected if a suitable electronic 
transition frequency happens to be very close to wjv- 



Conclusion. We have suggested a new experimental 
scheme to excite the nuclear 229 fTh - 229m Th transition 
in Th + ions and to accurately determine its frequency. 
The scheme relies on an electronic bridge process that 
is driven by two incident laser photons whose sum fre- 
quency is resonant with the nuclear transition frequency. 
Using our previous calculations of the electronic energy 
level structure of Th + in the experimentally relevant en- 
ergy range [7[ , we have estimated the probability of the 
investigated two-photon electron bridge process. Assum- 
ing that the nuclear transition energy is close to 7.6 eV, 
we find that a nuclear excitation probability in the range 
of 10 s _1 can be obtained with only moderate laser power 
and bandwidth requirements. 



This work was supported by Australian Research 
Council. The work of S.G.P. was supported in part by 
the Russian Foundation for Basic Research under Grant 
No. 08-02-00460-a. The work of E.P. and C.T. was sup- 
ported in part by DFG through QUEST. 



[1] E. Peik and C. Tamm, Europhys. Lett. 61, 181 (2003). 
[2] V. V. Flambaum, Phys. Rev. Lett. 97, 092502 (2006). 
[3] L. A. Kroger and C. W. Reich, Nucl. Phys. A 259, 29 
(1976). 

[4] C. W. Reich and R. G. Helmer, Phys. Rev. Lett. 64, 271 
(1990). 

[5] B. R. Beck, J. A. Becker, P. Beiersdorfer, G. V. Brown, 
K. J. Moody, J. B. Wilhelmy, F. S. Porter, C. A. Kil- 
bourne, and R. L. Kelley, Phys. Rev. Lett. 98, 142501 
(2007). 

[6] E. V. Tkalya, Phys. Uspekhi 46, 315 (2003). 
[7] S. G. Porsev and V. V. Flambaum, Phys. Rev. A 81, 
042516 (2010). 



032504 (2010). 

[11] 1. 1. Sobelman, Atomic Spectra And Radiative Transitions 

(Springer- Verlag, Berlin, Heidelberg, New York, 1979). 
[12] V. A. Dzuba, V. V. Flambaum, and M. G. Kozlov, Phys. 

Rev. A 54, 3948 (1996). 
[13] V. F. Brattsev, G. B. Deyneka, and I. I. Tupitsyn, Bull. 

Acad. Sci. USSR, Phys. Ser. (Engl. Transl.) 41, 173 

(1977). 

[14] M. G. Kozlov, S. G. Porsev, and V. V. Flambaum, J. 

Phys. B 29, 689 (1996). 
[15] V. A. Dzuba, M. G. Kozlov, S. G. Porsev, and V. V. 
Flambaum, Zh. Eksp. Teor. Fiz. 114, 1636 (1998), 
[Sov. Phys.-JETP 87 885, (1998)]. 



[8] URLhttp://www. lac. u-psud.fr/Database/Tab-energy/Th4WJtuEy.TBferib7cMi.Rtmkhnson, and P. Shorer, Phys. Rev. A 
[9] K. Zimmermann et al. (to be published). 26, 19 (1982). 

[10] S. G. Porsev and V. V. Flambaum, Phys. Rev. A 81, 



